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ABSTRACT 

The nonlinear decay of a disturbance was calculated by an iterative solution to the 
Navier- Stokes equations. The general expressions obtained for successive approxima- 
tions were suitable for machine computations, within the limitations of the computer ca- 
pacity. For the initial condition, a three-dimensional cosine distribution with two har- 
monic terms was assumed. The nonlinear interaction of these harmonic terms then 
produced new harmonics. As the iteration process proceeded, a large number of har- 
monics, at wave numbers higher than, lower than, or the same as the original ones 
were generated. This process appeared to be similar to that occurring in the early 
stages of the development of turbulence by flow through a grid. 
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NONLINEAR DECAY OF A DISTURBANCE IN AN UNBOUNDED VISCOUS FLUID 

by Robert G. Deissler 
Lewis Research Center 

SUMMARY 

The nonlinear decay of a disturbance was calculated by an iterative solution to the 
Navier- Stokes equations. The general expressions obtained for successive approxima- 
tions were suitable for machine computations, within the limitations of the computer ca- 
pacity. For the initial condition, a three-dimensional cosine distribution with two har- 
monic terms was assumed. The nonlinear interaction of these harmonic terms then 
produced new harmonics. As the iteration process proceeded, a large number of har- 
monics, at wave numbers higher than, lower than, or the same as the original ones 
were generated. This process appeared to be similar to that occurring in the early 
stages of the development of turbulence by flow through a grid. 


INTRODUCTION 

A fundamental problem in fluid mechanics is to predict what happens to a disturb- 
ance in a viscous fluid. Besides being of interest in itself, the problem is important as 
a first step in understanding the dynamics of turbulence. In particular, it should give an 
indication of the early stages of the development of turbulence from a regular disturb- 
ance, such as that produced by the flow of a fluid through a grid. 

Taylor and Green (ref. 1), starting with an initial three-dimensional sine-cosine 
velocity distribution, obtained successive approximations for the velocity as a function of 
position and time by substituting known approximations into the Navier- Stokes equations. 
The results, which showed that small eddies can be generated from larger ones, were 
limited to small times and to the first few approximations for velocities. Jain (ref. 2) 
obtained similar results by using the initial conditions from reference 1 in a numerical 
finite-difference solution of the Navier- Stokes equations. 

In this analysis, an initial three-dimensional cosine velocity distribution, which is 
more convenient to work with than that of Taylor and Green, is used to obtain iterative 
solutions of the Navier-Stokes equations. Of particular importance is the fact that these 



initial conditions can be conveniently used with the compact Cartesian-tensor form of the 
Navier- Stokes equations. Analytical expressions are first obtained'for the first three 
approximations. In each approximation, the linear terms of the Navier- Stokes equations 
are considered as unknown, and the nonlinear terms as known from the preceding ap- 
proximation. General expressions are then obtained which give the terms in an approxi- 
mation of any order in terms of those in the preceding approximation. By using these 
expressions, a computer can be programmed to generate successively higher-order ap- 
proximations from lower-order ones. It will later be seen that, in this process, new 
harmonics or eddy sizes are generated from existing ones. As the order of the approxi- 
mation increases, the Reynolds number for which the solution is applicable increases. 
Accordingly, the upper Reynolds number for which the solution is applicable will be lim- 
ited only by the computer capacity. Thus, the scope of the solution that can be feasibly 
obtained will be increased over that which could be obtained by carrying out the itera- 
tions algebraically by hand. The basic equations for the analysis will be considered in 
the next section. 


BASIC EQUATIONS AND INITIAL CONDITIONS 


The following is considered: given the initial velocity distribution in an unbounded 
viscous fluid, predict the motion at later times. The equations of motion for a viscous 
fluid with constant properties are, in dimensionless form, 


3Uj 8^u. 
dt dx^dx^ 


9p a < U i U k> 
3x. 3x k 
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where 
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The quantity U? is a velocity component, x? is a space coordinate, Xq is a character- 
istic length, t* is the time, p is the density, and v is the kinematic viscosity. (Note: 
stars on dimensional quantities are omitted on corresponding dimensionless quantities. 
Symbols are defined in the appendix. ) The subscripts in equations (1) and (2) can take on 
the values 1, 2, or 3, and a repeated subscript in a term indicates a summation. To ob- 
tain an explicit equation for the pressure, the divergence of equation (1) is applied to the 
continuity equation (2) to get 


3 2 p = _ a2 ( u f u k> (3| 

3 x i 3x je ax i 3x k 

In the remainder of the analysis, it will be convenient to use the set of equations (1) 
and (3), rather than (1) and (2). 

The expression assumed for the initial disturbance is, in dimensionless form, 

u^ = a? cos q • x + b? cos r" • x (4) 

where 


a? = ^ a?* 
1 v 1 


b° = 3“ b 9* 

1 i, 1 


q = x Q q * 


r 


= 


— * 
r 


0 * 0 * 

a^ and b^ are initial velocity amplitudes or Fourier coefficients of the disturbance, 

and q and r are initial wave number vectors. The quantities q • x and r ■ x are 
dot products (i. e. , q • x = q^x^ + q 2 x 2 + qgXg). 111 order to satisfy the continuity condi- 
tion (2), the a's, b's, q's, and r's in equation (4) must be related by 


3 



and 


(5) 


a ?«i = 0 




> 


Equation (4) could be generalized to include any number of similar harmonic terms. 
It represents the simplest expression that can give nonzero pressure and inertia terms 
when substituted into equations (1) and (3). It can be readily verified that, if only one 
term of equation (4) were used, the continuity conditions (5) would require the pressure 
and inertia terms (terms on right side of eq. (1)) to be zero, so that nonlinear effects 
would be absent. Because the present study is mainly concerned with nonlinear effects, 
two harmonic terms are included in equation (4). 


ANALYSIS FOR FIRST THREE APPROXIMATIONS 

As a first approximation for u^ at any time 

u. = a^ cos q • x + b^ cos r • x 


( 6 ) 


where a. and b^ are functions of time, the forms for which will be determined from 
the equations of motion. The form of equation (6) has been chosen so that it will reduce 
to equation (4) for t = 0. 

Substituting equation (6) into the right side of (3) yields 


d 2 ( u l u k ) 

9 V x k 


= cos ^ ~ ? ) • * - ^j? b j? r k a k cos ^ + 


• x 


(7) 


For obtaining equation (7) and subsequent equations, the following identities are used: 


sin A sin B = — (cos(A - B) - cos(A + B)]] 
2 

(8) 

sin A cos B = — [sin(A - B) + sin(A + B)1 
2 

(9) 
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cos A cos B = i [cos(A - B) + cos(A + B)] 
2 


( 10 ) 


In order to be consistent with equations (3) and (7), it is assumed that 

p = a cos(q - r ) • x + b cos(q + r ) • x (11) 


Substituting equations (7), (11), and (5) into (3) and equating coefficients of like cosine 
terms yield 


q 2 - 2q £ r £ + r 2 
q 2 + 2q £ r £ + r 2 


( 12 ) 


(13) 


where q and r are the magnitudes of the initial wave number vectors q and r (i. e. , 
q 2 = q^). Substituting equations (6) and (8) to (13) into the right side of (1) gives 


a P a < u i u k> 

ax i ax k 


_J_J (q. - r t ) + i (a^b* - b.r^) 


q 


2( lf r £ + r 


sin(<f - r) • x 


^(! b (? r k a k 1 

(q. + r t ) + i (syi^ + b^a*) 


q 2 + 2q^r £ + r 2 


sin(cf + r) ■ x (14) 


In order to be consistent with equations (1), (4), and (14), it is assumed that, for the 
second approximation for u. , 

u i = a i cos q ' * + b i cos r" • x + c^ sin(q - r") • x + d^ sin(q + r) ■ x (15) 


From the continuity condition (2) and equation (15), 



b.r. = 0 


c j( q i - r.) = 0 


djCqj + r i ) = 0 


To solve equation (1) by iteration, it is assumed that the linear viscous term 
3 2 u •/ 3x^3x^, as well as 3u^ /at, is unknown at each stage of the iteration. Considering 

3 2 u./ ^ x k^ x k as un ' :aiown should give a more accurate answer than would be obtained by 
considering it as known from the preceding approximation, as has been done in previous 
work. It will be shown that the present scheme gives a series of negative exponentials 
in time, rather than a power series, as was obtained in reference 1. Substituting equa- 
tions (15) and (14) into the left and right sides of equation (1), respectively, gives, after 
equating coefficients of like sine and cosine terms, 


Saj 

at 




= 0 


Using the initial condition (4) yields 


0 -qn 

a i = a i e 4 


(16) 


Similarly, 


b, . b?e- 2 ' 


(17) 


The use of equations (16) and (17) in the right side of equation (1) (i. e. , eq. (14)) 
then gives 


0 Cj 

at 


1 + (q 2 - 2q k r k + r 2 )c i 


0 „0 


-r t u HafV°-b° r 4 


q 2 - 2q jf r jf + r 2 


a -(q 2 +r 2 )t 


(18) 
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or, again by use of initial condition (4), 


c i = 


4q>r 


kk 


2a k r k b i q jg 
q 2 - 2q £ r £ + r 2 


(q. - r.) + afb^ - b°a°r A 


X 


e -(q 2 +r 2 )t 


-(q 2 -2q k r k +r 2 )t 


(19) 


Similarly, 


d i = 


4 % r 


kk 


2 4 r k b K 

q 2 + 2q £ r £ + r 2 


<«i + r t> + a ? b « q J + b ? a « r t 


X 


e -(q 2 +r 2 )t 


-(q 2 +2q k r k +r 2 )t 


( 20 ) 


The third approximation is obtained in similar fashion by substituting the second approx- 
imation into the right sides of equations (1) and (3). 

The first, second, and third approximations for u. may be summarized as follows: 

First approximation: 


u i = a f e_C1 * cos q • x + b?e 


0 -rn 


cos r • x 


( 21 ) 


Second approximation: 


7 



2 2 

u i = a i* e ~^ * cos Q • x + b?e _r * cos r • x 


4r, q 


k^k 


2a k r k b jt q l / x 0.0 .0 0 

(q.-rp + ^b^-b.a^r^ 


q 2 - 2q £ r^ + r 2 


.(q 2 + r 2 )t . -<q 2 -2V fc+ r 2 )t 


sin(q - F) • x 


4 qir r 


kk 


2a k r k b £ q £ (q . + r . } + a 0 b 0 q ^ + b?a 0 r£ 


q 2 + 2q £ r £ + r 2 


X 


" -(q 2 +r 2 )t r2 >‘ 

6 c 


sin(q + r ) • x (22) 


Third approximation: 



(Cont. on next page) 
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(Cont. on next page) 
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4q t r 


k A k 


******** < q i + r i) + a K q k + b ? a k r k 


q 2 + 2q k r k + r 2 


e -(q 2 + r 2 )t . e -(q 2 +2q k r k+ r 2 )t 


sin(q + r ) • x 


a?r 


~ fCkqk ■ < 2 «i - r i> - ; (“Mk - c ? a k r k 


4q 2 - 4q k r k + r 2 


x < 


2q - 4q k r k L 


,(-2q 2 -r 2 )t _ e -(4q 2 -4q k r k+ r 2 )t 


2q " 2q k r k L 


(- 2q 2 +2q k r k - r 2 )t - (4q 2 - 4q k r k +r 2 )t 

6 - G 


cos(2q - r ) • x 




4q 2 + 4q k r k + r 2 


< 2 «i + r i> - \ (^ d K - d i’ a k r k) 


x 


2q + 4q k r k L 


e (-2q 2 -r 2 )t . e -(4q 2 +4q k r k+ r 2 )t 


(Cont. on next page) 
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1 


2q 2 + 2 q k r k 


^(-2q 2 -2q k r k -r 2 )t 


-(4q 2 +4q k r k +r 2 )t 


cos(2q + r ) • x 


c l r A 

q2 - 4q k r k + 4r 2 


(9 i - 2r i> - ; - b ? c k r t) 




; 4 Vk + 2r ' 


1 (-q 2 -2r 2 )t _ -(q 2 -4q k r k+ 4r 2 )t 


2 Vk - 2r ‘ 


(-q 2 +2q k r k -2r 2 )t -(q 2 -4q k r k +4r 2 )t 

e - e 




>cos(q - 2r ) - x 


d S r ‘- 0 


— - (4i + 2r i> - ; ( d ? b K + b i d k r k) 


q 2 + 4q k r k + 4r 2 


.Vk + 2r ‘ 


e (-q 2 -2r 2 )t _ e -(q 2 +4q k r k+ 4r 2 )t 


2q k r k + 2r L 


(-q 2 -2q k r k -2r 2 )t ^-(q 2 +4q k r k +4r 2 )t 


-T'v 


J] 


cos(q + 2?) • x 


(Cone, on next page) 
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2r i 0 ,0 / 0,0 ,0 0\ 

- — c $ T ll d k T k + { c i d k + d i c k) i 


3 (-2q 2 -2r 2 )t _ -4r 2 t 


2 2 
-q + r 


(-2q 2 -2q k r k -2r 2 )t _ 4 2 fc (-2q 2 +2q k r k ~2r 2 )t _ 4 ^ 

e - e e - e 


-2q 2 - 2q k r k + 2r 2 


-2q 2 + 2q k r k + 2r 2 


sin 2r • x 


2c ^i 0 0 

— c A d k q k 
Lq 2 


/ 0,0 ,0 o\ 

■ V c i d k + d i c k) 


q k 


1 (-2q 2 -2r 2 )t _ -4q 2 t 


2 2 
q - r 


(-2q 2 -2q kV 2r 2 )t _ ^ (-2<, 2 + 2 Vk -2r 2 )t _ ^ 


2q 2 - 2q k r k - 2r 2 


2q 2 + 2q k r k - 2r 2 


sin 2q ■ x 


(23) 


where a? and b? are given by the initial condition (4), and 


c° = 


1 

2a k r k b t<lj> 

4 Vk 

q 2 - 2q £ r £ + r 2 


<4i - r i> + a ? b f«f - b i a l r l 


(24) 


<f 


1 

*Vk 


2 4 r A 

q 2 + 2q J r 1 + r 2 


<«i + r i> 


a l b ' °iy 


.0 0 

+ b i a 4 r j(l 


(25) 


It is evident that the number and complexity of terms in the expression for u. in- 
crease rapidly as the order of the approximation increases. This happens because each 
term in a particular approximation is multiplied by every other one when it is substituted 
into the nonlinear terms of equations (1) and (3); that is, each sine or cosine term inter- 
acts with every other one to produce more terms. If, on the other hand, the nonlinear 
terms in equations (1) and (3) were neglected (low Reynolds number), the first approxi- 
mation (eq. (21)) would be an exact solution of (1) and (3). 
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GENERAL EXPRESSIONS FOR NONLINEAR PRODUCTION OF HARMONICS 


Because of the proliferation of new harmonic terms as the order of the approxima- 
tion increases, it appears impracticable to generate approximations much higher than 
the third analytically. Rather, it would seem better to obtain general expressions for 
the new terms produced by the interaction of pairs of harmonic terms in a particular ap- 
proximation. It should then be possible to program a computer to generate successive 
approximations from preceding ones. These general expressions for the nonlinear inter- 
action of pairs of harmonic terms will be considered next. 

Examination of the terms in the first three approximations for u^ (eqs. (21) to (23)) 
shows that they are of the general form 


u. 


2 K J cos 7 • X + A? - sin k ■ x) 

K 


(26) 


where 


and 



(27) 


The superscripts s and c appear in coefficients of sine and cosine terms, respectively. 
The indices k , q, and r take on as many values as are required to represent the terms 
of the particular approximation. 

When equation (26) is substituted into the nonlinear terms of equations (1) and (3), 
each term of equation (26) is multiplied by (or interacts with) every other term. Thus, 
sine - sine, sine - cosine, and cosine - cosine products must be considered. Consider, 
first, the terms generated by the interaction of two terms 


and 


A. c — cos n • x 
i, n 

> 


A, C - 

i, m 


cos in • x J 


(29) 


13 


Thus, the sum 


AP - cos n ■ x + A? — cos m ■ x 


(30) 


where 



(31) 

(32) 


is substituted for u. in the right sides of equations (3) and (1). The manipulations are 
similar to those used for obtaining the first and second approximations in the last sec- 
tion. Because it is necessary to obtain the additive contribution to u^, which results 
from the interaction of the two harmonic terms in equation (30), u. is replaced on the 

Of* * 

left side of equation (1) by u which represents that contribution. For the initial 

i,n, m’ 

condition on u c? -S u c, _S — is set equal to zero when t = 0. (Initial condition (4) will 
i, n, nr i, n, m 

be satisfied later when the totality of terms contributing to u. is considered. ) The ex- 

pression for u c? -S _ is, then, 
l, n, m 


> r- 


U Cj i * = 
n, m 



2 m f a i,n ) q\ a k, m, r ( . - ) + a c - n ap _ - a c _ m. a£ - , 

2 9 1 1 z, n, q kk, m, r 1 , m, r k k, n, c 


2 2 
n - 2n k m k + m 




n 2 - 2rum, + m 2 ~ b~ - bS. 

k k m, r n, q 


sin(n - in) • x 



( 

q, r 


2m a. c m - 


5 , r (n . + m.) + a c _ n a, c - + a c _ m a, c _ 

2 1 1 i,n,q k k,m,r 1, m, r k k, n, q 


n + 2n k m k + m‘ 


e -fe, r +b n-, q ) t . e -( n2+ \ m k +m2 ) t _ 

K c sin(n + m) ■ x 


n 2 + 2n,m, + m 2 - bL - b£. 
k k m. “ 


n, q 


(33) 
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Similarly, the interaction of terms like 


and 


A sin n • x 
i, n 


A _ sin m ■ x 
i, m 


( 34 ) 


where 


and 


Z -bi t 

af _ e n ’ q 
1 , n,q 


( 35 ) 


results in 


« = V af — 

1, m / y 1, m, : 


-bl t 

m, r 


( 36 ) 


u s ’ i _ 



2m af n af — 
j g 1, n, q k k, m, r 


2 2 
n - 2n^m^. + m 


(n. - m.) + a s n a s - af — m a 
x i i, n, q k k, m, r i, m, r k k, n, q 


q, r 


-(b® +bl It 
n 7 q m, r 


-(n 2 - 2 n k m k+ m 2 )t' 


n 2 - 2n, m. + m 2 - bf - bf* 

k k n, q m, r -J 


sin(n - m) • x 



2m af _ n af — 

1 £ , n, q k k ? m, r 

2 2 
n + 2n k m^ + m 


(n. + m.) - af n. af — - af — nr af — 

x i' i, n, q k k, m, r l, m, r k k, n, q 




■( b ff,q* b a,r)'. -( A 2 W m2 )' 


n 2 + 2n^m k + m 2 - b^ - bf* 


sin(n + fn) ♦ x 


( 37 ) 


n, q m, r -J 
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Finally, the interaction of terms like 




A — sin n • x 
i, n 


and 


A. cos m • x 
1 , m 


J 


(38) 


where A s — and A c are defined by equations (34) and (30), results in 
i, n i, m 


u s ^ 
i, n, m 



2m a® _ n, a, c - 


h* (n . - m.) + a s _ n, a, c - - a c _ m a® _ 

? 2 i r i, n, q k k, m, r 1 , m, r k k, n, i 


o 2 

n - 2n k m k + m 


q,r 




n 2 - 2n, m, + m 2 - bl - 

k k n, q m, r J 


cos(n - m) • x 



2m a® _ n a, c _ 


— Q ^ — (n. + m.) - a® — n a c — - a c _ nia® _ 

2 2 i 1 l, n, q k k, m, r 1 , m, r K k, n, q 


2 2 
n + 2n k m k + m 


<b r 


-(bl +b c L V -(n 2 +2n rm +m 2 )t 
5 \ n, q m, r/ \ k k / 

n 2 + 2n. m. + m 2 - bl - b^. 

k k n, q m, r 


cos(n + m) • x 


(39) 


The successive approximations to the velocity components Uj may now be calcu- 
lated. If the v th approximation is known, equations (33), (37), and (39) give contribu- 
tions to the v + 1 st approximation, which are due to the presence (or interaction) of 
pairs of harmonic terms in the approximation. Thus, if the v^ 1 approximation for 
Uj is given by 


Id 



(40) 


= . . . + A. c _ cos n • x + Af sin if • x 
1 1 , n l, n 


+ . 


+ Ap cos in • x + A s _ sin in • x + . 
x, m i, m 


then the 


v + 1 


st 


approximation is 


,(v+l) _ 


,c ? c 


+ U. + 


+ u 


S, s 


+ . 


. + U S ’-S + 


+ u s? ^ + . . . + aPe ^ ^ cos cf • x + b?e~ r ^ cos r • x (41) 

i, m,n 1 1 

where the uj L. _Js are given by the equations (33), (37), and (39). 

The last two terms in equation (41) are the first approximation (eq. (21)) and are in- 
cluded in equation (41) because, whereas every other term is regenerated at each stage 
in the iteration process, the terms of the first approximation are not so regenerated. 

That is the case because of the assumed initial condition u Cj -P = 0 when t = 0. The 

i, n, m 

cc 0 0 _► i 

same result could be obtained by letting u. = a. cos q • x + b. cos r • x at t = 0 

J l, n, m i i 

when the two interacting terms include the first approximation. However, it is simpler 
to use equation (33) throughout as it stands and to include the first approximation in equa- 
tion (41). 

In obtaining equation (41) from (40), the presence (or interaction) of each pair of 

harmonic terms in equation (40) produces a term u^ ^ — in (41). Thus, each term in 

i, m, n ’ 

equation (40) interacts once with every other term. Because of the continuity condition 

(2) (see also eq. (5)) no contributions are obtained by the interaction of a term with it- 

th st 

self. Thus, if there are s terms in the v cn approximation, the v + I s approximation 
will have (s - 1) + (s - 2) + . . . + 1 = (s/2)(s - 1) contributions, like u^ L. plus the 
first approximation. 

With the velocity components u. known to some degree of approximation the 

2 1 

space-averaged values of uf can be calculated. Because u. is given by equation (26), 
squaring, integrating over a cycle, and using the orthogonality property of sines and co- 
sines result in 
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RESULTS AND DISCUSSION 

If the velocity distribution in a fluid is given by a sum of harmonic terms as in equa- 
tion (26), then the presence (or nonlinear interaction) of each pair of terms produces new 
harmonics given by equations (29) to (39). Each pair of terms then present interacts to 
produce still more harmonic terms or eddy sizes. It is evident that there will result a 
very large number of eddy sizes, or a situation resembling turbulent flow. It should be 
pointed out that, because of wave number vectors if - m and n + m, which occur, for 
instance, in equation (33), the generated harmonics can have either higher or lower wave 
numbers than the original ones. Moreover, when some of the generated harmonics in- 
teract with the original ones, new harmonics can be generated with the same wave num- 
bers as the original ones. 

One might ask how a random turbulence could develop from a regular disturbance 
produced, for example, by flow through a grid with equal spacings. It has been pointed 
out earlier that because of continuity, harmonics of the same wave number and/or inten- 
sity vectors will not interact nonlinearly. However, if there are small random disturb- 
ances in the flow upstream of the grid or small imperfections in the grid, the eddies 
produced by the grid will not be identical; therefore, they can interact to produce new 
harmonics. Because of the random nature of the small upstream distrubances, or grid 
imperfections, the resulting flow will be random or turbulent. 

Figure 1 shows calculated kinetic energy per unit mass for a disturbance plotted as 
a function of dimensionless time. First, second, third, and fourth approximations are 
shown. The initial condition for the disturbance is given by equation (4), where the 

values of a?, b? q. and r, are shown in the figure. The values of the initial vectors 

1111 0 0 
were chosen to satisfy continuity and to give nonzero interaction terms (a. q. = b r. = 0, 

0 0 iiii 

a i r i ^ anc * ^i^i ^ 0). Disturbances having three different initial intensities are com- 
pared. The first approximation (eq. (21)) gives the linear decay that occurs when the 
nonlinear or inertia terms in the Navier-Stokes equation are neglected. Comparison of 
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Kinetic energy per unit mass, 


12xl0 3 



0 . 005 . 010 . 015 . 020 . 025 . 030 . 035 . 040 

Dimensionless time, t 


(a) Weaker disturbance Initial dimensionless Fourier coefficients: 
aj = 20(1, -1, -2,) ; hip- 80Cl f -2, -1). 


45xl0 3 


24xl0 3 



0 . 005 . 010 . 015 . 020 . 025 . 030 . 035 . 040 

Dimensionless time, t 

(b) Moderate disturbance. Initial dimensionless Fourier coefficients: 
af=30(l, -1, -2); b® - 120(1, -2, -1). 

Figure 1. -Time decay of mean kinetic energy per unit mass. 
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Dimensionless time, t 

(c) Stronger disturbance. Initial dimensionless Fourier coefficients: 
aP- 40(1. -1, -2); bf = 160(1, -2, -1). 

Components of initial wave number vectors: qj = (4, 2, 1); rj = 2(4, 1, 2). 


the nonlinear with the linear solutions shows that the nonlinear effects cause the disturb- 
ances to decay more slowly. It is noted that a similar result is obtained for random 

turbulent flow. Strongly nonlinear turbulence decays approximately as t~*, whereas the 

-5/2 

decay law for linear turbulence is closer to t ' . 

Comparison of the various approximations indicates that at large times, higher 
order approximations are required for the iteration process to converge (as indicated by 
the agreement of the third and fourth approximations) than are required at smaller times. 
Larger values of the initial intensity parameters indicate convergence only at small 
times for the first four approximations. On the other hand, smaller values of the initial 
intensity show convergence even at very large times. The poorer convergence of a given 
approximation at larger intensities is, of course, due to the stronger inertia effects that 
occur at high Reynolds numbers; therefore, higher approximations are required to take 
them into account. 

To show how the energy is distributed among the various directional components, 
and Ug 

case shown corresponds to that in figure 1(b). Comparison of the curves in figure 1(b) 
with those in figure 2 shows that the third and fourth approximations in the curves in fig- 
ure 2 converge for a smaller range of dimensionless times than do those of figure 1. 
Evidently, a given order approximation at a time t is less accurate for the individual 
components than it is for the total energy. Comparison of the first (linear approxima- 
~2 

tion for Ug 

with time because the different components decay viscously at different rates, depending 
on the initial conditions. 


j Uj and for u u^ shows that the ratios can either increase or decrease 


the ratios 


2/2 

U 2/l 


A 


are plotted against dimensionless time in figure 2. The 


Approximation 


First (linear) 

Third 



Figure Z - Variation of mean square velocity component 
ratios with time. Initial dimensionless wave number 
vectors: qj *(4,2,1); r ? = 2(4, 1,2). Initial dimension- 
less Fourier coefficients a^ = 30(1, -1, -2); 
bp = 120(1, -2, -1). Moderate disturbance. 
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The fourth (nonlinear) approximations show the same general trends as the first approxi- 
mations, but the nonlinear effects tend to bring the ratios closer to one (this is also true 

for Ug/ Ug); that is, the nonlinear effects tend to make the disturbance more isotropic. 

To show how various harmonics contribute to the kinetic energy of the disturbance, 
kinetic energy spectra were calculated from equation (42) and plotted against dimension- 
less wave number k in figure 3. The cases shown here are the same as those in fig- 
ure 1. Results are given for several dimensionless times for which the iteration process 
has converged (see fig. 1). For t = 0, there are, of course, only two spikes or har- 
monics in the spectrum because of the assumed initial condition (eq. (4)). For the later 
times shown, particularly for the stronger disturbances, a large number of other har- 
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(a) Weaker disturbance. Initial time, t = 0. 
Initial dimensionless Fourier coefficients: 

20(1, -1, -2); bP- 80(1, -2, -1). 
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(b) Weaker disturbance. Time, t = 0.02. Initial dimensionless Fourier coef- 
icients: a^ = 20(1. -1, -2); bf- 80(1, -2, -1). 
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(c) Moderate disturbance. Initial time, t - 0. 
Initial dimensionless Fourier coefficients: 
aJ-30U f -1, -2); bf - 120(1, -2, -1). 


(d) Moderate disturbance. Time t = 0.02. Initial dimensionless Fourier coef- 
ficients; a? - 30(1, -1, -2); bp = 120(1, -2, -1). 


Figure 3. - Spectrum of mean kinetic energy per unit mass. Components of initial wave number vectors: qj - (4, 2, 1); r { - 2(4, 1, 2). 
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(e) Stronger disturbance. Initial time, t = 0. 
Initial dimensionless Fourier coefficients; 
ap 40(1, -1, -2); bP= 160(1, -2, -1). 
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(f) Stronger disturbance. Dimensionless time, 1 = 0.01. Initial dimension- 
less Fourier coefficients; ap = 40(l r -1, -2); bp= 160(1, -2, -1). 


Figure 3. - Concluded. 


monies have developed because of the interaction between contributions at various wave 
numbers. Although some of the generated harmonics are strong, many of them are too 
weak to have an appreciable effect on the disturbance, so that their representations lie 
close to the horizontal axis. By increasing the strength of the disturbance, more of the 
harmonics will have appreciable strengths. However, the fourth approximation will then 
not be sufficient to obtain convergence of the iterative process. Unfortunately, the com- 
puter used for the calculations did not have enough storage capacity to carry out the ap- 
proximations higher than the fourth. 

For comparison with the nonlinear results just discussed, the spikes for the linear 
solution are shown dashed in figures 3(b), (d), and (f). In this case, no new harmonics 
are generated. The original harmonics simply decay independently. Comparison of the 
solid and dashed lines shows that the harmonics at the original wave numbers are greatly 
affected by nonlinear interactions; that is, many of the eddies, generated or destroyed 
by nonlinear effects, have the same wave numbers as those in the original disturbance. 

It might be mentioned that, in order for that to occur, it is necessary to go to at least 
the third approximation because the second approximation produces only contributions 
that are at higher or lower wave numbers than the original ones (eq. (22)). 

Although the harmonics generated by nonlinear effects can occur at either high or 
low wave numbers, most of the significant generated harmonics in the calculated results, 
particularly those for the weaker disturbances, tend to lie in the low wave number re- 
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gion. The difficulty of generating significant harmonics at the higher wave numbers ap- 
pears to be related to the fact that viscous effects in the high wave number region tend to 
prevent their formation. As the disturbance strength increases, however, the tendency 
to form significant harmonics in the higher wave number region appears to increase. 


CONCLUDING REMARKS 

An initial three-dimensional disturbance consisting of two harmonics can decay non- 
linearly and generate new harmonics at wave numbers greater than, smaller than, or the 
same as, those of the original harmonics. If either the wave number vectors or the am- 
plitude vectors of the two initial harmonics have the same direction, however, the conti- 
nuity condition requires that there be no interaction. In that case or if nonlinear effects 
are ignored by neglecting the inertia terms in the Navier- Stokes equation, the original 
harmonics will simply decay independently. 

In general, nonlinear effects cause the disturbance to decay more slowly. At low 
Reynolds numbers there is, of course, no significant difference between the linear and 
nonlinear solutions. 

The ratios of the directional components of the mean square velocities changed with 
time. This was caused by differences in the rates of viscous (linear) decay of the vari- 
ous directional components and by transfer of energy between the components by non- 
linear action. Whereas the viscous effects could change the component ratios in either 
direction (depending on the initial conditions) the nonlinear effects always tended to bring 
the ratios closer to one or to make the disturbance more isotropic. 

For a given approximation, the iterative solution for the nonlinear case was most 
accurate at small times and/or small Reynolds numbers of the disturbance. As either 
time or the Reynolds number increased, the order of the approximation that was required 
to obtain convergence of the iterative process increased. 

The generation of a range of harmonics in the early stages of the development of 
turbulence by flow through a grid appears to be similar to the production of harmonics in 
the present analysis. Random turbulence can evidently develop by flow through a grid 
because of small random disturbances in the flow upstream of the grid or small imper- 
fections in the grid. Random harmonics are then generated as a result of the interaction 
of the slightly different intensities and wave numbers of the disturbances in the flow from 
the grid. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 9, 1968, 

129-01-05-20-22. 
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APPENDIX - SYMBOLS 


A c A s A c 

' a, k’ i, K 1 , n’ 

A c A s A s _ 

1, m i, n i ? m 

A c *_, A s *_ 

1 , K 1 , K 


a, b 


dimensionless Fourier coefficients for velocity components. 


Fourier coefficients for velocity component 


dimensionless Fourier coefficients for pressure fluctuation, 
X 0 a */P^ 2 , x o b */ py2 (see eq. (11)) 


a i> b i> c i> d i 


ap — , a s _ , a c _ , 

i,K,q’ i,K,r’ i, n, q 

CSS 

a. _ , a. _ , a. _ 

1, m, r ? i, n, q’ i, m, r 

a 0 b° 

i 5 u i 


0 * , 0 * 
a. . b. 

i ’ i 


a*, b* 

a*, b* 


dimensionless Fourier coefficients for velocity component, 
x o a i*A> vrA etc - (see eqs. (6) and (15)) 

constants given by eqs. (27), (28), (31), (32), (35), and (36) 


initial dimensionless Fourier coefficients for velocity com- 
0* / 0* / 

ponents, x Q a i jv and x-b. / v (see eq. (4)) 
initial Fourier coefficients for velocity components 

Fourier coefficients for velocity component 

Fourier coefficients for pressure fluctuation 


b- ,bl , bS. , 

K, q K, r’ n, q 

bl , bl , bl 

m, r n, q’ m, r 
m, n 


m F n i 


m*, n* 


constants given by eqs. (27), (28), (31), (32), (35), and (36) 


magnitudes of dimensionless wave number vectors 
components of dimensionless wave number vectors, 
x 0 n i* 

components of wave number vectors 


m, if 
P 


q i’ r i 


dimensionless wave number vectors 
dimensionless pressure, x^p * / pv^ 
pressure 

components of initial dimensionless wave number vectors, 
x 0 q i*’ x 0 r i* 
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components of initial wave number vectors 


q* r* 
n ’ x 


q, r 

t 

. * 


u. 

1 

* 

u i 


u. 


.C, c 


,u s ’« - 


1 , n, nr 1 , n, m 1 , n, m 


x i 

X* 

1 


x 

X r 


K 

V 


P 

Superscripts: 

c 

s 

* 


initial dimensionless wave number vectors 

9 

dimensionless time, ^t*/ x Q 
time 

component of dimensionless velocity, x^u:* / v 
component of velocity 

interaction terms given by eqs. (33), (37), and 

dimensionless position coordinate, x^/Xq 
position coordinate 
dimensionless position vector 
characteristic length 

dimensionless wave number vector, x^/T* 
wave number vector 
kinematic viscosity 
density 

cosine terms 
sine term 

dimensional quantities 
averaged quantities 
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